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Abstract 

A modified method for calculating the pion vacuum susceptibility from an effective quark-quark 
interaction model is derived. Within this approach it is shown that the pion vacuum susceptibility 
is free of ultraviolet divergence and is much smaller than the previous estimations. 
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The vacuum condensate is known to play the essential role for characterizing the non- 
perturbative aspect of the QCD vacuum. In fact, the gauge invariant vacuum condensates 
such as the chiral quark condensate, the two gluon condensate and the mixed quark gluon 
condensate are well known examples studied extensively so far [1,2]. Meanwhile susceptibili- 
ties of vacuum are also important quantities of strong interaction physics. They directly enter 
in the determination of hadron properties in the QCD sum rule external field approach [3-5]. 
In particular, the strong and parity-violating pion-nucleon coupling depends crucially upon 
X 71 " , the 7r susceptibility [6]. The previous estimation for the value of pion susceptibility were 
obtained by a QCD sum rule three-point formalism [7] . It is the aim of this letter to con- 
sider 7r vacuum susceptibility in the framework of the global color symmetry model(GCM). 
Though GCM[8] violate local color SU(3)c gauge invariance and renormalizability, they pro- 
vide a successful description of various nonperturbative aspects of strong interaction physics 
and hadronic phenomena at low energies [8-22]. A few more recent review articles on this 
topic can be found in Refs. [23,24]. 

The GCM generating functional for massless quarks in Euclidean space can be written 
as [8] 

Z[fj, 77] = j VB d (x, y) eMS% V, B e (x, y)}) , (1) 

with 



S[fj,r,,B e (x,y)] = -Tr]n[p5(x-y) + ±A 9 B B (x,y) 



+ 



J Jd 4 xd 4 y 



B e (x,y)B e (y,x) g 

+ r](x)g(x,y; [B \)r](y) 



(2) 



2gW(x-y) 

where B e (x,y) is an auxiliary bilocal field introduced as in Ref. [8] and the matrices A e = 
D a F b C c is determined by Fierz transformation in Dirac, flavor and color space (more detail 
can be found in Ref. [8]); D ab u (x — y) denotes the effective gluon two-point green function. 
Because the form of D ab v {x — y) in the infrared region is unknown, one often treat the 
D ab u (x — y) as the GCM input parameter, which, as we will discuss later, is chosen to 
reproduce certain aspects of low energy hadronic physics. For convenience we will use the 
Feynman like gauge D ab v [x — y) — S^ u S ab D(x — y) from now on. 

The saddle-point of the action is defined as SSlfj, rj, B e (x, y)]/5B e (x, y) = and is 

r]=r]=0 

given by 

B e (x -y) = g 2 D{x - y)tr[\ e G {x - y)\, (3) 
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where Q stands for G[Bq] and the trace in Eq.(3) is to be taken in Dirac and color space, 
whereas the flavor trace has been separated out. 

We will calculate the vacuum condensates from the saddle-point expansion, that is, we 
will work at the mean field level. This is consistent with the large Nq limit in the quark 
fields for a given model gluon two-point function. In the mean field approximation, the field 
B e [x — y) is substituted by their vacuum Bq{x — y). Under this approximation, the dressed 
quark propagator G(x — y) = Qo(x — y) in GCM has the decomposition 

G-\p) = i 7 • p + n(p) = i 7 • P A(p 2 ) + B(p 2 ) (4) 

with the self-energy dressing of the quarks II (p) through the definition: 

"I 



n(p) = ^A 6 B 6 (p) = J d 4 xe ipx 



t 1 -p[A(p 2 )-l] + B(p 2 ), (5) 



where the self energy functions A{p 2 ) and B{p 2 ) are determined by the rainbow Dyson- 
Schwinger equation[8] 

In order to get the numerical solution of A(p 2 ) and B(p 2 ), one often use model forms for 
gluon two-point function as input in Eqs.(6,7). As a typical example, we choose g 2 D(q 2 ) = 
4:7r 2 <i^+A with d — and three sets of different parameters for A and A, which are adjusted 
to reproduce the pion decay constant in the chiral limit f n =87MeV and the chiral low energy 
coefficients[14,25] . 

Here we want to stress that the B(p 2 ) in Eqs.(6,7) has two qualitatively distinct solutions. 
The "Nambu-Goldstone" solution, for which 

B(p 2 ) ^ 0, (8) 

describes a phase in which: 1) chiral symmetry is dynamically broken. Because one has a 
nonzero quark mass function; and 2) the dressed quarks are confined, because the propaga- 
tor described by these functions does not have a Lehmann representation [9]. In "Nambu- 
Goldstone" phase, the vacuum configuration B°(x — y) in GCM (at the mean field approx- 
imation) can be regarded as a good approximation to the "exact" vacuum in QCD. The 



alternative "Wigner" solution, for which 

B(p 2 ) = 0, (9) 

describes a phase in which chiral symmetry is not broken and the dressed- quarks are not 
confined. 

With these two "phase" characterized by qualitatively different quark propagator, the 
GCM can be used to calculate the vacuum condensates and susceptibilities as we will show 
below. 

In the external field of QCD sum rule two-point method, one often encounters the quark 
propagator in the presence of the J r (y) = q(y)Tq(y) current (r stands for the appropriate 
combination of Dirac, flavor and color matrices). 

s% T {x) = (o\T[q c a(^(o)]b)jr = s% v > PT {x) + s% v > NP {xi (10) 

where S^ ,PT (x) is the quark propagator coupled perturbatively to the current and ST ' NP (x) 
is the nonperturbative quark propagator in the presence of the external current J r (one 
should note that the external current should be taken to be J r (f>r, where 0r is the value of 
the external field. In what follows, to simplify the notation, we will take the 0r=l, which 
does not affect our results). The vacuum susceptibility x r m the QCD sum rule two-point 
external field treatment can be defined as [7] 

S% T ' NP {x) = (0| : &(x)$(0) : |0) jr = —T^S^ H(x)(0\ : 9(0)9(0) : |0>, (11) 

where the phenomenological function H(x) represents the nonlocality of the two quark non- 
local condensate. Note that H(0) = 1. 

The presence of external field implies that S^ r (x) is evaluated with an additional term 
AL = — J r • r added to the usual QCD Lagrangian. In three point method of QCD sum 
rule, if one takes only a linear external field approximation, the S^ /3 r,NP (x) in Euclidean 
space is given by [7] 

S C a f' NP (x) = J dty0| : q c a {xmy)Tq\y)q^) : |0>. (12) 
Using Eq.(ll), (12) and note that if (0) = 1 we have 

- ^5 cc T^ X r (0| : 9(0)9(0) : |0) = / d 4 y(0\ : q :(0)q e (y)T q e (y)q^0) : |0). (13) 
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Multiplying Eq.(13) by T Pa 5 cc >, we get 

16tt 2 



r 

X a 



J d 4 y(0\ : g c (0)rg c (0)g e (y)rg e (y) : |0>, 



(14) 



tr 7 (rr) 

with a=-(27r) 2 (0| : <?(0)g(0) : |0). Eq.(14) shows that the vacuum susceptibilities originates 
from the nonlocal four quark condensate [7]. Therefore, in order to calculate the vacuum 
susceptibilities, one must know how to define and calculate the vacuum condensates in 
advance. 

Let us first recall the definition of vacuum condensate in QCD sum rule. In QCD sum rule, 
one often postulate that quark propagators are modified by the long-range confinement part 
of the QCD; but the modification is soft in this sense that at short distance the difference 
between exact and perturbative propagators vanishes. 

To formalize this statement, one can write the "exact" propagator G(x) as a vacuum 
expectation of a T-product of fields in the "exact" vacuum |0) 



G«(x,i/) = (0|T[ ft (x)«,-(y)]|0>. 



According to the Wick theorem, one can write the T-product as the sum 



(15) 



T[qi(x)qj 



Qi(x)Qj(y) + : Qi{x)qj{y) ■ 



(16) 



of the "pairing" and the "normal" product. The "pairing" is just the expectation value of 
the T-product over the perturbative vacuum |0) 



qi {x)q 3 {y) = (0\T[q t (x)q 3 (y)}\0) = G^ rt (x,y). 



(17) 



i.e., the perturbative propagator. By this definition, we have the following two quark vacuum 
condensate (0| : q{x)q(y) : |0): 



(0| : qi (x) qj (y) : |0> = (0|T[fc (x)qj(y)] |0> - (0\T[q t (x) qj (y)}\0) 
= (-) \G ji (y,x)-G p j r t (y,x)] = UVufax) = (-) 



(18) 



A Jq-(y-x) 



(27T) 



G{q 2 ) - G pert {q 2 ) 
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where E(x, y) = G(x,y) — G pert (x ) y). 

Similarly, we have the four quark vacuum condensate: 



i : q(x)A^q(x)q(y)A^q(y) : |6> 

(0|T [q(x)A^q(x)q(y)A%(y)} |0) - (0|T [q(x)A^q(x)q(y)A^q(y)} |0> 



(0| : q{x)A^q{x)q{^A^ q {y) : |0) - (0\: q(x)A^q(x)q(y)A^q(y) :\0) 
(0| : q(x)A^q(x)q(y)A^q(y) : |0> - (0| : q(x)A^q(x) q(y)A™q(y) : |0> 



- {tr[E(y, x)A (1) S(a:, y)A (2) ] - ir[S(x, x)A (1) ]tr[S(y, y)A (2) ]} , 



(19) 



here the A^ stands for an operator in Dirac and color space. 

Based on the above statement, in order to calculate quark vacuum condensates, one must 
know not only the "exact" but also the "perturbative" quark propagator in advance. The 
calculation of "exact" quark propagator in GCM(at the mean field approximation) has been 
given by Eqs.(6,7) with Eq.(8), the key question left now is how to treat consistently the 
"perturbative" quark propagator in GCM. Because the form of perturbative quark propaga- 
tor in the low momentum region is unknown(one only know the large momentum behavior of 
perturbative quark propagator), it is practically very difficulty to provide a G pert for whole 
momentum region in QCD. However, in phenomenological application, one may proceed by 
making some assumptions regarding the "perturbative" quark propagator. 

As a first "guess" for a perturbative quark propagator, one would simply take the bare 
quark propagator(A(p 2 )=l, B(p 2 )=0) or the perturbative quark propagator at one loop ap- 
proximation. However, numerical calculation for the mixed quark gluon condensate in GCM 
shows that one can not get a reasonable numerical value for the mixed quark gluon conden- 
sate from this simple "guess" . In order to obtain the reasonable result for the mixed quark 
gluon condensate, the authors in Ref.[26] identified the "perturbative" vacuum with the 
"Wigner" vacuum(no chiral symmetry breaking and the dressed-quarks are not confined), 
in contrast to the "Nambu-Goldstone" vacuum corresponding to dynamical symmetry break- 
ing and the dressed-quarks are confined. 

In "Wigner" phase(the B'(p 2 ) = 0), the Dyson-Schwinger equation(6,7) reduces to: 



where the solution A'(p 2 ) in Eq.(20) denotes the quark propagator propagating in the 
"Wigner" vacuum(It also sums infinitely many effective one-gluon exchange diagrams, as 
does the solution in "Nambu-Goldstone" phase). If one use the vacuum configuration B'® of 
the "Wigner" mode(at the mean field approximation) to model the "perturbative" vacuum 
in QCD (As it did in Ref.[26]), the "perturbative" quark propagator in GCM can be written 
as G pert (q 2 ) = -jrPzf 2 = —i / J-qC(q 2 ). It should be noted that the "perturbative" propagator 




(20) 
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here is only a notation(In order to be consisted with the one used in QCD sum rule) which 
denotes the dressed quark propagator in "Wigner" phase. Numerical studies also show that 
for q 2 3> A 2 QCD) one has 



(21) 



A 2 (q 2 )q 2 + B 2 (q 2 ) 

i.e., at short distance the difference between exact and perturbative propagators tends to 
zero(see Eq.(27) below). This result is just what one expected in advance. It can be seen 
from the following figures (fig.l-fig.2). 



Fig.1 . Comparison between A(s) and A'(s) 
for gluon propagator g 2 D(q 2 )=4ji 2 d^. 2 /(q 4 +A) 



Fig.2.B(s) for gluon propagator g 2 D(q 2 )=47t 2 dA 2 /(q 4 +A) 



3.5 



2.5 



1.5 



A(s) setl (A=0.1GeV 4 ,?i=1.77GeV) 
" " A'(s) setl 

- A(s) set2 (A=0.01GeV 4 ,^=1.33GeV) 
1 - 1 A'(s) set2 

A(s) set3 (A=0.0001GeV 4 ,X=0.95GeV) 

— A'(s) set3 





s (GeV 2 ) s (GeV 2 ) 

Once the "exact" and "perturbative" quark propagator in our model are determined, one 
can calculate the (0| : qq : |0), (0| : qA^qqA^q : |0) vacuum condensate at the mean field 
level. In particular we obtain the two quark condensate (0| : qq : |0) in the chiral limit; 

3 (V- , B(s) 



(0|:gg:|0) = (-)tr[S(x,0)]| 



^-Hj™ — J o dss-^ + B2{s) . 



(22) 
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It should be noted that Eq.(22) is appropriate only when we use an effective gluon prop- 
agator; i.e., g 2 D(q 2 ), falls off like 1/g 4 or more quickly (for example using the form of 
g 2 D(q 2 )=A7r 2 d-^ r ^). In this case the model does not need renormalization and the chiral 
quark condensate defined in Eq.(22) is free of UV divergence (There are many such models, 
the instanton liquid model among them). However, if we take into account the asymptotic 
behavior of the gluon propagator (for example using the form of Eqs. (25,26) below), the 
renormalization is necessary, and the naive formula Eq.(22) for the condensate has to be 
multiplied by a renormalization constant(more detail can be found in Ref.(12) and (16)). 

Returning to the calculation of vacuum susceptibilities, in the case of tensor current (r = 
<T^ U ), from Eqs. (14), (19) and the fact tr^a^a^) = 48, we obtain the tensor vacuum 
susceptibility x T a: 



T 

X a 



(48tt 2 



d A p 



B{p 2 



A 2 (p 2 )p 2 + B 2 (p 2 ) 



3 lim 

JO 



sds 



B(s) 



A 2 (s)s + B 2 (s) 



(23) 



The calculation of the tensor vacuum susceptibility in the framework of GCM has been given 
in Ref.[27]. 

Similarly, in the case of pseudoscalar current, we have the pion susceptibility x^a: 




A{p 2 



1 2 



A 2 (p 2 )p 2 + B 2 (p 2 ) 
A(s) 



C( P 2 



C(s) 



p 



s + 



B(p 2 



1 2 



A 2 (p 2 )p 2 + B 2 (p 2 ) 

2' 



A 2 (s)s + B 2 (s) 



(24) 



A 2 (s)s + B 2 (s) 

which is UV convergent. It is reasonable due to the vacuum condensates reflecting 
the IR behavior of QCD. However, if we do not adopt an adequate subtraction mecha- 
nism^. e., C(p 2 )=0) to calculate the n vacuum susceptibility, Eq.(24) would be UV diver- 
gent. In addition, it should be noted that the calculation of nonlocal four quark condensate 
(0| : q c (0)Tq c (0)q e (y)Tq e (y) : |0) here is quite different from the previous estimation [7]. In 
Ref. [7], the authors did not adopt subtraction mechanism to calculate nonlocal four-quark 
condensate, instead they used vacuum saturation hypothesis and retained only the contri- 
bution of the scalar condensates(more details can be found in Ref. [7]). 

From Eq.(22) and (24), one can calculate the two-quark condensate and the pion vacuum 
susceptibility. In order to show the importance of the proper choice of the perturbative 
quark propagator, we use separately the "Wigner" solution of Eq.(20) and the bare quark 
propagator for the perturbative propagator to calculate the pion vacuum susceptibility x^a 



and x l7Ta m Eq.(24). In Table. I, the results for the two quark vacuum condensate and ir 
vacuum susceptibilities for three different parameters sets of the model gluon two-point 
function are displayed. 



TABLE I: The two quark condensate and ir vacuum susceptibility for g 2 D(q 2 )=Air 2 d^ 



A[GeV 4 } 


X[GeV] 


-(gg)5 (MeV) 


X w a (GeV 2 ) 


X ,7T a (GeV 2 ) 


HT 1 


1.77 


291 


0.109 


6.402 


1(T 2 


1.33 


250 


0.089 


4.603 


10- 4 


0.95 


217 


0.071 


3.583 



Tabled show that our two quark condensate is compatible with the "standard" value of 
— (250Me\/) 3 in QCD sum rule, whereas the ir vacuum susceptibility x" a is much smaller 
than the value x? a — (1-7 — 3.0) GeV 2 obtained within a phenomenological approach [7]. In 
addition, the huge difference between the y^a and x' Wa clearly shows the importance of the 
proper choice of the perturbative quark propagator. Due to the subtraction of bare quark 
propagator(or the perturbative quark propagator at one loop) is not appropriate for the 
calculation of the mixed quark gluon condensate [26], we take the x* a as our prediction for 
the 7i vacuum susceptibility. 

As is shown in Ref.[26], in a correct treatment of vacuum condensate, the cutoff \x should 
be taken to be infinity. Because the calculation is numerical, so one have to use a very 
large but finite value of upper limit of integration. In our numerical calculation, the upper 
limit is chosen to be 500 GeV 2 , which is large enough to ensure that the calculated value 
is independent of the choice of the upper limit of integration. This result can be seen from 
Table II. 

To check if the small value of the pion susceptibility is due to a special form for the 
effective gluon propagator, other two popular forms of the gluon propagator[14,25] have 
been employed. They are model 1: 

^ V) = Ml(?) + g'Duvtf) = Sir^e^ + ^ ^ (25) 

\ QCD / 
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and model 2: 



^V) = <?D?ltf) + g 2 D uv tf) = A^d^— + v . (26) 



\ 1S -QCD 



TABLE II: The sensitivity of the vacuum condensates to the choice of the upper limit of integration 
H in Eqs.(22) and (24) for g 2 D(q 2 )=4TT 2 d^ K (A = lQ- 2 GeV 4 , X = 1.33GeV) 



H (GeV 2 ) (-)(0| : qq : |6)3 (MeV) x w a (GeV 2 ) 
1 179 0.083 
10 241 0.089 
100 249 0.089 
300 250 0.089 
500 250 0.089 



The term Dm(q 2 ), which dominates for small q 2 , simulates the infrared enhancement 
and confinement. The other term Duv{q 2 )i which dominates for large q 2 , is an asymp- 
totic ultraviolet (UV) tail match the known one-loop renormalization group result with 
d = [12/(33 - 2N f )] = 12/27, A QCD = 200 MeV. The model parameters A and A are 
adjusted to reproduce chiral low energy coefficients and the pion decay constant in the 
chiral limit as previously did. The calculated results are shown in Table III. 

TABLE III: x 71 a f° r Models 1 and 2 with three sets of different parameters, respectively. 



Model 1 


Model 2 


A(GeV 2 ) 


X(GeV) 


X n a(GeV 2 ) 


A(GeF 4 ) 


X(GeV) 


X n a(GeV 2 ) 


0.200 


1.65 


0.096 


io- 1 


1.83 


0.112 


0.020 


1.55 


0.066 


io- 4 


1.02 


0.069 


0.002 


1.45 


0.058 


io- 7 


0.83 


0.081 



From Table. Ill, one can see that the pion vacuum susceptibilities are not sensitive to the 
model gluon propagators and the small value in comparison to Ref. [7] is robust within this 
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model approach. Here we want to stress that Eq.(24) is convergent even if we take into 
account the asymptotic behavior of the gluon propagator. This is quite different from the 
calculation of chiral quark condensate(22). Taking into account the asymptotic behavior 
of the gluon propagator in Schwinger- Dyson equation, for q 2 — * oo, we have the proper 
asymptotic behavior for quark propagator 

^-^''''- (ffM-^ (27) 

which, numerically, will appear very close to zero(see Fig.l and Fig. 2) because it is powerlaw 
suppressed. However, for the chiral quark condensate in Eq.(22), we have 

<6| : qq : |0> oc jT sds ^^ oc {l W }^ , (28) 

which is logarithmic divergent (j m is the mass anomalous dimension). From Eq.(27), it is 
easy to check that Eq. (24) is free of the UV divergence. This means that using the forms of 
Eq.(25) and (26) do not cause any problems in the calculation of pion vacuum susceptibility. 

To summarize, in the present letter, a subtraction mechanism which is consistent with 
the definition of vacuum condensate (the vacuum expectation value of the normal produc- 
tion of quark operator) in QCD sum rule, is adopted to calculate pion vacuum suscepti- 
bility(employing Wick theorem and considering in addition to the usual Nambu-Goldstone 
phase a perturbative phase of GCM). Within this approach the pion vacuum susceptibility 
is free of UV divergence. This is reasonable and what is to be expected. The results of pion 
susceptibility we obtained are much smaller than the previous estimation, which can be used 
as a reference for further study. Finally, we want to stress that GCM is not renormalizable. 
Therefore, the scale at which a condensate is defined in our approach is a typical hadronic 
scale, which is implicitly determined by the model gluon propagator g 2 D(q 2 ) and the solu- 
tion of the rainbow DS equation(6,7). This situation is very similar to the determination 
of vacuum condensates in the instanton liquid model where the scale is set by the inverse 
instanton size [28,29]. 
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